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Abstract 

The data from solar neutrino experiments together with standard solar model 
predictions are used in order to derive the possible profile of the magnetic field 
inside the Sun, assuming the existence of a sizeable neutrino magnetic moment and 
the resonant spin flavour mechanism. The procedure is based on the relationship 
between resonance location and the energy dependent neutrino suppression, so that 
a large neutrino suppression at a given energy is taken to be connected to a large 
magnetic field in a given region of the Sun. In this way it is found that the solar 
field must undergo a very sharp increase by a factor of at least 6-7 over a distance 
no longer than 7 - 10% of the solar radius, decreasing gradually towards the surface. 
The range in which this sharp increase occurs is likely to be the bottom of the 
convective zone. There are also indications in favour of the downward slope being 
stronger at the start and more moderate on approaching the solar surface. Typical 
ranges for the magnetic moment are from a few times 10 -13 /i£ to its laboratory 
upper bounds while the mass square difference between neutrino flavours is A 2 m,2i — 
(0.6 - 1.9) x W- 8 eV 2 . 
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1. Introduction 



It is usually argued that the solar neutrinos which were first observed in 1968 [|TJ have 
not provided so far, as opposed to previous expectations, any particular knowledge of the 
physics inside the Sun. In fact the solar interior is opaque and detailed information about 
its structure has turned out possible to obtain, not via direct probes, but indirectly, via 
helioseismology [0. 

Neutrinos however hold the potential of real probes of the solar interior, since they 
can travel across the Sun practically unperturbed. A neutrino produced near the centre 
of the Sun typically has a collision probability of the order of 10~ 6 — 10 -7 while travelling 
through the Sun, so the energy spread of a neutrino beam while traversing solar matter 
is left unaltered. If neutrinos are endowed with magnetic moment, an idea originally 
proposed by Cisneros |§ and later revived by Voloshin, Vysotsky and Okun J|, they can 
be converted via spin flip while travelling through the solar magnetic field into weakly 
sterile particles. Due to their neutrality the direction of propagation is undeflected. 

One of the original motivations for this hypothesis was the apparent anticorrelation 
of neutrino flux in the Homestake experiment with sunspot activity ||. The data seemed 
to suggest that the more and the larger the sunspots, the fewer neutrinos were detected, 
indicating a more efficient spin flip conversion at times of more intense solar activity. 
Although this is an attractive hypothesis, several statistical analyses fail to confirm 

this possibility and the data from all experiments [1C] - [14] appear to be consistent with 
a constant rate Jl5[ . 

It is now widely acknowledged that there exists a solar neutrino deficit which, in very 
general terms, may be stated as too few neutrinos being detected fT(| - |14| as compared to 



theoretical predictions |Td| - [£1]]. The possibility of a solution to this discrepancy within 
the minimal standard electroweak model appears highly unlikely. Hopes for solutions 
based on wrong interpretations of the experiments |TD[ - [|TJJ] or different input solar 
physics have successively been abandoned ||22j| . More and more the right solution seems 
to lie on non-standard neutrino properties such as oscillations (both MSW 



or vacuum 



the resonant spin flavour precession mechanism (RSFP) |25| or flavour changing 
neutral current interactions [2(J. The RSFP mechanism combines the original idea of 



Cisneros P| with the interesting possibility of the conversion mechanism being a resonant 
one, like matter oscillations |23j and is totally consistent, as will be seen, with the absence 
of the above mentioned anticorrelation. 

The present paper deals with RSFP type solutions to the solar neutrino problem. 
This mechanism requires a smaller magnetic field than the simple spin flip and one of its 
features is the fact that neutrinos of different energies have their resonances in different 
locations within the Sun, an essential point in the present analysis. 

In the past few years of solar neutrino physics, it has become increasingly clear that the 
deficit is energy dependent p7] - |33j , in the sense that weakly active neutrinos of different 



energies are suppressed differently. Since the resonance range is the zone where most if 
not all the conversion occurs and a stronger magnetic field leads to a greater suppression, 
an energy range where neutrinos are largely suppressed is a clear indication of a large 
magnetic field in a certain range of the Sun's interior. Conversely, a small suppression in 
a given energy range indicates a small magnetic field in the region where the corresponding 
resonance is located. This principle is the basic idea of the present work and in this way 
a 'tomography' of the solar magnetic field may be obtained. Hence the neutrino energy 
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suppression spectrum is expected to be directly connected to the profile of the magnetic 
field along the Sun's radius, provided the neutrinos have a sizeable magnetic moment. It 
should be stressed that a sensible solution for the field profile is reasonably constrained, as 
it must account for a strong suppression for the intermediate energy neutrinos, a moderate 
one for the energetic S B ones and no suppression for the pp ones. This result together 
with the requirement that all neutrinos should have a resonance inside the Sun implies, 
as an essential feature, that the average magnetic field should undergo a discontinuity or 
at least a very fast variation in space somewhere in the interior of the Sun. 

Recent progress in helioseismology seems to confirm the fact that the fluid character 
of the convective zone is interrupted near its bottom at approximately 0.71-Rs. Below 
this depth a rigid body like structure appears. This is thus the likeliest region for the fast 
variation to arise. From the analysis of this paper one can predict that, starting from the 
solar surface where the field should not exceed a few kG in the sunspots, one expects an 
increase along the convective zone towards a maximum near its bottom. This increase 
appears to become more intense in the deeper layers of the convective zone. Proceeding 
inwards in the radial direction, a sudden drop by an order of magnitude or more over the 
following 7 — 10%i2s is expected. Below this range the field is either negligible or may 
rise very moderately. This expectation is consistent with a neutrino magnetic moment in 
the range from a few times 10~ 13 /i£ up to its laboratory bound of 10.8 x 1CT 10 /ib @- 

The paper is divided as follows: in section 2 a brief review of the experimental sit- 
uation - [14] is presented and its implications for the survival probabilities of low, 



intermediate and high energy neutrinos are derived using the predictions from several 
standard solar models |U| - In section 3 the neutrino propagation equation through 
solar matter is solved analytically and generalized to include the neutrino energy [33] and 
production zone distributions and non-adiabatic effects in the form of the Landau 



Zener approximation [^| . The expression of the survival probability is next evaluated. In 
section 4 the necessary general profile which is required in order to satisfy experimental 
data and the results from solar models is established. It is argued that the high suppres- 
sion found for the intermediate energy neutrinos and the absence of suppression for the 
low energy ones along with the fact that their energies are relatively close is an indication 
of the sudden field variation referred to above. An example of a toy magnetic field is 
obtained from approximate general conditions. In section 5 an outlook is presented. Ma- 
jorana neutrinos with transition magnetic moments [y e — > v^) will be used throughout, 
the results for Dirac neutrinos being practically identical ||37|| . 



2. Energy Dependent Solar Neutrino Suppression 



The data from the four solar neutrino experiments |10[ - |14j are quoted in table I and 
the corresponding predictions from six solar models [16|] - |2T| including the contribution 
from each neutrino component are given in table II. All models include heavy element 
diffusion except for TCL []I(J and TCCCD | I7| . It is now generally acknowledged that a 
'standard' solar model (SSM) should include diffusion, owing to the fact that such models 
give a remarkably good agreement with data from helioseismology 0. 

The absence of the intermediate energy neutrinos, consisting principally of the 1 Be 
line at E = 0.8QMeV and the CNO continuum, has been realized several years ago [2D] 
from the comparison of the Homestake and Kamiokande data. It is considered by some as 
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the 'true' solar neutrino problem in the sense that it is independent from normalization to 
any solar model, either standard or non-standard f38|| . It appears as a natural consequence 
of the luminosity constraint plj , |32[| (Xq — 1.367 x W^Wcm' 2 ) 

L Q = Y,(®-< E »>k)fk (k = pp,pep, 7 Be,CN0 8 B) (1) 

k 

with Q = 26.73MeV (total energy released in each neutrino pair production) and the 
equations [[3^] 

S Ga = E^i/i (t = PP,pepjBe,CNO, 8 B) (2) 

i 

Sci = E^cijj (i= 7 Be,CNO, 8 B) (3) 
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f pep = 0.0023/ pp (4) 

where Sea, Sci denote the Gallium and Chlorine data. 

One has in fact from this system of four equations, upon elimination of the pp flux 



f pp , using the average value Sa a — 73.8 ± 7.7SNU and the nuclear cross sections [|32jj 
a Ga,ii &ci,ji the following intermediate energy neutrino flux (in units cm _2 s _1 ) 

f Be = 1.04 x 10 4 /b - 2.88 x 10 10 (5) 
fcNO = -8.46 x 10 3 f B + 2.22 x 10 10 . (6) 



Inserting f B from SuperKamiokande, |T1] the total flux from these neutrinos is negative: 

f Be = -3.42 x lO^cm-V 1 (7) 
f CNO = 1.56 x KPcm-V 1 . (8) 



Better fits were done by the authors of |32|,[3Ij] who obtained 

fse+CNO < 0.7 x lO^m- 2 ^ 1 (3a) (9) 
f Be+ CNO = (-2.5 ±1.1) x lflWV 1 . (10) 

These authors used the former Kamiokande flux data which were higher than Super- 
Kamiokande. From equations (5) and (6) it is seen that the total flux fse+CNO decreases 
with decreasing f B , so that the results (9), (10) should be further aggravated in the non- 
physical direction. Hence the probability that neutrinos are standard is no greater than 
1%, while, if the luminosity constraint is dropped, it may increase to 4% pi]] . So inter- 
mediate energy neutrinos appear in practice to be completely suppressed. 

As far as high energy ( 8 B) neutrinos are concerned and denoting by Rci the ratio 
data/SSM prediction one may write 

Rc^R^Pi + R^Ph- (11) 

Here Rq x is the fraction of intermediate (high) energy neutrinos in the Chlorine experi- 
ment as theoretically predicted and Pi(h) is the fraction of intermediate (high) energy v eL 
produced in the Sun that are detected on Earth. Using Pj = (99%CL) and the models 
listed in table I one gets for P H the range 

0.35 <P H < 0.63 (12) 
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with the smaller value corresponding to BP95 |2TJand the larger to TCCCD |T7J . 

The same argument may now be applied to the low energy (pp) neutrinos. With 
obvious notation, one has for the Gallium result 



R Ga = R L Ga P L + R^Pj + R% a P H (13) 

Again, taking Pj = 0, one obtains for the two extreme cases 

P L = 0.981 (BP95)[21] (14) 
P L = 0.966 (TCCCD) (15) 

which means that the low energy (pp) neutrinos are practically unsuppressed. Thus any 
solution to the solar neutrino problem based on non-standard neutrino properties must 
be consistent with the following general conclusions: 

1. No suppression for low energy (pp) neutrinos. 

2. Strong suppression for intermediate energy (J Be, CNO) neutrinos. 

3. Moderate suppression for high energy ( 8 B>) neutrinos. 

The requirements Pl > 0.96, Pj < 0.04, 0.35 < P# < 0.63 will be used throughout the 
paper. 



3. Survival Amplitudes and Probabilities 

The starting point of this section is the Schrodinger-like evolution equation for Majo- 



rana neutrinos in the Sun |}7[] 



( ^ ^ 


M 








»m 1 



(16) 



with the Hamiltonian H given by 




11 - I ' "..n A 2 m 21 , (G=)lN e • ' ' ' ' ) 

Here A 2 rri2i is the mass square difference between mass eigenstates, N e is the electron 
density taken as exponentially decreasing with radius in the region of interest |35| , GN e = 
2.11xl0~ u exp(—r /0.09Rs)eV , B is the solar magnetic field and the remaining notation is 
obvious. A vanishing vacuum mixing angle is assumed throughout and the mixing angle 
from matter and magnetic moment/field is [57] 



,,2 R2 

6 ^B* + \[^ 2 N e -^ + ((I^N e -^y + A^B^Y 

The neutrino propagation through solar matter is assumed to be adiabatic except for a 
short range |37j in the neighbourhood of a critical point - the resonance - where the mixing 
angle (18) is maximal so that the density verifies 

5G ; ^ = ^ (19) 



3V2 e 2E 
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Across the resonance the neutrinos may shift from one eigenvalue to another, this effect 
being parametrized by a jump amplitude. So a neutrino produced as a weak eigenstate 
v eL at Ti which is a superposition of two mass eigenstates v\, v 2 



VliXi) \ _( C6i S0i\(\ 



(20) 



will become at location n, after having passed through a resonance [3{J, the following 
combination 



"lin) \_ ( Xiexp i J^ 1 H Dll dr X 2 exp i J£ H D22 dr \ ( v x {r i} \ 
vii?\) J I -X 2 exp i J r ri H Dll dr \\exp i J r n H D22 dr I I u 2 (r. A '' 



Here Ai denotes the amplitude for adiabatic propagation, A2 is the jump amplitude across 
the resonance and Hu n , H^ 22 are the 11 and 22 elements of the diagonalized Hamiltonian. 
Denoting by 81 the diagonalizing angle at r 1; it is straightforward to derive the survival 
amplitude, that is, the amplitude that a neutrino which is produced as v eL remains as v eL 

Amp(is eL (ri) -> v eL (n)) = ce^Xiexp i H Dll dr - c dx s 6i \ 2 exp i H D22 dr 

+s 01 cg i X 2 exp i H Dll dr + se^e^lexp i H D22 dr. 

(22) 

For the conversion amplitude one just performs in (22) the replacements cq 1 — > s$ 17 —sg 1 — > 
cq 1 . In a simplified calculation where all neutrinos are considered to be emmited at 
the same point and with the same energy, the survival probability is simply the square 
modulus of eq. (22) where all fast oscillating terms are set equal to zero. This is the well 



known Parke's formula 39 



Psurv — ^ + ~ -Pc)c26» l C 2 6» 1 (23) 

where the quantity Pq is the jump probability, Pc = IA2I 2 , which will be taken in the 
Landau Zener approximation |36|, f|0| []. 



Pc = expj-TT^^OMRs). (24) 

2E 



The coefficient of this exponent is the adiabaticity parameter, 



lc = 2 -^0.09R s . (25) 



A 2 r»2i 
2E 



A more accurate calculation must include the production and energy range of the 
neutrinos. Denoting by x(E) and 0(r«) the production amplitudes per unit energy and 
unit length respectively, the survival amplitude at the edge of the Sun (ri = Rs) becomes 
the sum of the following four terms 



1 It was previously shown that the Landau Zener approximation which assumes a linearly decreasing 
density in the vicinity of the critical point works rather well in the Sun. See |4(| 
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= / x(^) c 0i / \i4>(ri)c ei exp i / H Dll drdndE (26) 

J E m J a Jri 

[Em rb rRs 

■ I x( E ) c 6i I \ 2 (j)(ri)se,exp i / H D22 drdridE (27) 

J E m Jo, JTi 



rE M rb rR s 

+ / x{E)se 1 j \24>(ri)c ei exp i / H Dll drdridE (28) 

J E m J a Jri 

rE M rb rRs 

+ / x( E )s 01 / \\4>{ri)s0 i exp i / H D22 drdridE. (29) 

J Em Jo, Jri 



while the conversion amplitude 

rE M rb r R s 



r^M ro rns 

/ x( E ) s 0i / \i<p{ri)c ei exp i / H Dll drdridE (30) 

J E m Jo Jri 



Rs 

X{E)s 6l / \ 2 4>(r i )so i exp % \ H D22 drdridE (31) 



1M 

/ x(E)c dl / \* 2 (j)(ri)ce t exp i / H Dll drdridE (32) 

/ x(E)c 01 / Ai0(rj)s ei exp z / H D , 22 drdridE. (33) 

J _E m , J a Jri 



In these expressions the first two integrals in each term extend over the energy and 
production range for each neutrino component (pp, 7 Be, 8 B, ...) as given in refs. |35[, [^1|] 
respectively (model with He and metal diffusion). 

Turning now to the probabilities, it should first be observed that nuclei produce neu- 
trinos in a disordered manner both in space-time and in energy. Each neutrino therefore 
acquires a random phase endowing the terms (26) - (29) and (30) - (33) which therefore 
add incoherently, as opposed to (22). The random neutrino emission thus makes the 
amplitudes completely incoherent, so one adds probabilities. The generalization of (23) 
including energy and production range distribution is thus 



rE M rb 

Psurv = g(E) f(n 

J Em. J & 



C2 e C2 e dridE. (34) 



For the conversion probability one just performs in (34) the replacement (+) — > (— ). 
The functions g(E), f(ri), representing respectively the probability density of neutrino 
production per unit energy and per unit length, will be taken from the numerical tables 
in refs. [J35| , pi ]. An example of a survival probability showing the characteristic 'bath 
tub' shape is shown in fig. 1. 

Equations (34), (24), (25) constitute the basis of the following analysis in section 4. 
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4. Magnetic Field Profiles 



4.1 General Features 

Apart from the knowledge that neutrinos may provide us if they have a sizeable mag- 
netic moment, very little is known about the solar magnetic field. From the observation 
of the interstellar field, a lower bound may be fixed at IO/iG, whereas an upper bound of 
(2 — 5) x 10 3 C7 near the surface may be considered as associated to the field in sunspots. 
There is however no indication as to the extension of these structures and its associated 
field in depth significantly below the surface layers. 

Several 'plausible' profiles have been proposed and investigated in the litterature pll - 



|44| but these analyses which are now several years old did not take into account an energy 
dependent solar neutrino suppression. It should first be emphasized that the resonance 
condition (19) fixes the location of the critical density for each neutrino according to its 
own energy |37| : 

A 2.11xl(T 11 eV A 
x res = Omiog 3 ^ A2 — (35) 

2E 

where x is the fraction of the solar radius x = r/Rs- 

The range around the critical density where the Hamiltonian eigenvalue splitting, 



cti — Ol2 



(36) 



is dominated by the term 4fi 2 B 2 is where most if not all the conversion (i. e. neutrino 
suppression) occurs. 

An efficient conversion means that the jump probability (24) across the resonance 
from one Hamiltonian eigenvalue to the other is small, so each neutrino remains in its 
own eigenvalue. The process is therefore adiabatic and jc (25) is large. Most of the 
neutrinos will gradually change their flavour and the space rate of this change is maximal 
at the vicinity of the critical point. Conversely, a sizeable jump probability means that 
a significant fraction of neutrinos will change from one Hamiltonian eigenvalue to the 
other, thereby staying in the same weak eigenstate and the process is non-adiabatic. In 
all resonant processes the resonance is of course located between the production zone and 
the surface of the Sun, so the mixing angle (18) at the production point is close to zero 
and at the edge of the Sun is close to n/2 [|3"7|| . In the rough approximation (23) with no 
energy range nor production zone distribution, this means C20i — 1 and C2e 1 — —1 so that 

Psurv^Pc (37) 

A large jump probability opposes conversion, thus favouring survival. In such case 
most neutrinos, while travelling across the resonance, will leave to the other Hamiltonian 
eigenvalue which corresponds, past the resonance, mainly to the initial weak eigenstate. 
So neutrino survival is related to strongly adiabatic processes and conversion to the lack of 
adiabaticity. Recalling the analysis of section 2, these arguments show that pp neutrinos, 
for which P^>0.96, must correspond to highly non-adiabatic resonances, 

7c<0.01 (38) 
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while intermediate energy neutrinos which are highly suppressed (P/<0.04) must corre- 
spond to 

7c>l- (39) 
Requiring condition (38) to hold for all pp neutrinos whose energies lie in the interval 

0.236MeV<E pp <0A20MeV (40) 

and taking intermediate energy neutrinos satisfying (39) to consist totally of the 7 Be line 

at 

Er Be = 0.86MeV (41) 
one has respectively from (38), (39) using (25) 

0.09R s 2fi 2 B 2 (x? Be ) > (42) 

0.09R s 2fi 2 B 2 (x pp ) < 0.01— ^i. (43) 

The solar magnetic field should satisfy these two conditions at the corresponding critical 
densities. Using the values for E7 Be = 0.86MeV and E pPmax = 0A2MeV, it is seen from 
(35) that the spacing between these critical densities is 

^X7 Be , pp = (0.065 - 0.116) (44) 

implying a sharp rise of the magnetic field by at least a factor of (6 - 7) over a relatively 
short radial distance. Therefore it is seen from (42), (43) that the field is both bound 
from below and from above at points whose mutual distance is Ax7 Bepp : 

5.22 x 10~ 3 r— , . 

B(x7 Be ) > v/A 2 m 21 (45) 



B(x pp ) < 



0.483 



/ \ E pp 



A 2 m 2 i , . 

(46) 



Here f is the neutrino magnetic moment in Bohr magnetons, \i v = ffiB, energies are in 
eV and B in Gauss. No information is provided from this analysis as to the location of 
this particular topology of the solar magnetic field. Helioseismology indicates that such 
a behaviour can only appear around the bottom of the convective zone at approximately 
0.71Rs- This constrains, as seen from (35), the parameter A 2 m 2 i to lie within less than 
an order of magnitude from 10~ 8 eV 2 . 

As far as 8 B are concerned, their spread in energy is much larger than the previous 
ones [ 35| 

0.814MeV < Es B < IhMeV (47) 

so some of their resonances overlap with the one from 1 Be. Given the survival probability 
expected for these neutrinos as derived in section 2 

0.35 <P H < 0.63, (48) 

one may naively expect from (37), throughout most of their critical density range, 7c to 
be situated in the interval 

70^(0.15- 0.30). (49) 
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However, due to the large energy spread here, no definite general criterion can be estab- 
lished for the field distribution along the critical density range. For a short interval at the 
beginning of this range owing to the overlap with the 7 Be resonance, the field should of 
course be constrained as discussed, but for the remainder, the actual profile is essentially 
unconstrained, except that the field should be substantially smaller on average. Only the 
total probability (12) can be used as an indicator, together with the requirement that the 
surface field should not exceed a few kG. It may remain quite large over a considerable 
extent of the critical density range, so that an almost total conversion will occur for a 
substantial fraction of 8 B neutrinos with 7c far exceeding in this zone the interval (49). 
Consequently, in the remaining part, the field should be much smaller with highly non- 
adiabatic transitions and 7c well below (49). Such is the case with one of the distributions 
to be examined below. So the requirement that 7c lies within the interval (49) for the 
resonances of S B neutrinos is an indicative one for distributions varying smoothly in the 
range x > xi Be - 



4.2 Numerical Examples 

The general characteristics of suitable magnetic fields expressed in inequalities (45), 
(46) will now be confronted with specific field distributions. Equation (34) for the survival 
probability will be used together with the requirements Pj < 0.04, Pi > 0.96 and 0.35 < 
Ph < 0.63 as discussed in section 2. 

The starting point is a field spread only over the convective zone proposed years ago 
Hfl. It changes from zero at 0.7Rs to 10 5 G at 0.75i?s, then remains constant until O.8-R5 



and falls linearly to 2 x 10 3 G at the surface of the Sun. For simplicity, a linear increase 
from 0.7 Rs will also be assumed: 



B(x) = 0<a:<0.7 

B(x) = 2 x 10 6 (x - 0.7) 0.7<z<0.75 

B(x) = 10 5 0.75<s<0.8 

B(x) = 10 5 - 4.9 x 10 5 (x- 0.8) 0.8<x<l. 



It is straightforward to realize that this field satisfies conditions (45), (46). Requiring 
in fact (45), (46) to hold simultaneously for the 7 Be and all pp resonances, i. e. up to 
E PPmax = 0A2MeV, the intensity of the field in which pp resonances are located is at most 
1/7 of the field at the 7 Be resonance. Since the distance between the critical densities 
corresponding to E pPmax and Ei Be is 0.065i?s, (see eq. (35)), this leads, using (50), to 

X7 Be - 0.065<0.707 (51) 

(i. e. all pp resonances should be located deeper than 0.707i?s) and 

X7 Be > 0.700, (52) 

this last condition to ensure a non-vanishing magnetic field at the 7 Be resonance. In- 
equalities (51), (52) imply using (28) 

8.0 x 10~ 9 e1/ 2 <A 2 m 2 i<1.8 x 10~ s eV 2 . (53) 

This double inequality thus reflects the compatibility between the pp and 7 Be solutions 
in (50) according to conditions (45), (46). No constraint from the S B flux probability 
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has yet been imposed. Using the lower bound of (53) in (45) with the corresponding 
critical density field B(x7 Be ) = 10 5 G to determine within this criterion the lowest neutrino 
magnetic moment compatible with pp and 7 Be solutions, one obtains 

H„ > 4.7 x 10- 12 fi B . (54) 

Using this value in eq. (34) for the survival probability with P c as in (24) one finds, 
however, a minor overlap in A 2 m 21 as regards pp and 7 Be neutrinos: 

8 B (4.0 - 5.95 x i(T 9 )eU 2 (9.0 x KT 8 - 1.25 x 10~ 7 )eU 2 

7 Be (1.2 x 10- 9 - 1.13 x 10- 8 )eU 2 (55) 

pp < 4.3 x 10- n eU 2 > 6.85 x 10- 9 eU 2 . 

This overlap [A 2 m 2i = (0.69 - 1.1) x 10" 8 eU 2 ] is due to the fact that (45), (46) were 
imposed in this example for all pp neutrinos which need not be the case. A small fraction 
of pp neutrinos may have their resonances in a very short range with a field slightly in 
excess of 1/7 of the field at the 7 Be resonance, provided a substantial fraction of pp 
resonances is located in a range where the field is much smaller. In terms of x pp this 
overlap range where (45), (46) are not satisfied is 

0.707 < x pp < 0.721 (56) 

and for smaller x pp the field (50) is of course below 1/7 of B( 7 Be). On the other hand 
the solution for S B neutrinos in (55) is clearly inconsistent with the other two. It can be 
turned however consistent by relaxing (45), (46) over a short range, hence relaxing (54). 
In this way the set of solutions allowed by (50) may be investigated. The result is 

3.1 x 10 _1 Vb<^<3-8 x 10' 12 fi B , 6.33 x 10~ 9 eU 2 <A 2 m 21 <6.56 x 10~ 9 eU 2 (57) 

respectively. In order to have an idea of the amount of maximum deviation from (45), 
(46) involved, take the 'worst' solution, namely 

fi v = 3.1 x lO" 12 ^ , A 2 m 2i = 6.33 x 10~ 9 eU 2 . (58) 

Condition (45) then gives 5(a;7 Be )>133970G' (actually B(x7 Be ) = 10 5 G) while (46) is 
satisfied for all neutrinos with energies E<0.3A2MeV . The field distribution (50), shown 
in fig. 2, is therefore a borderline case in terms of providing a solution to the solar neutrino 
problem. 

What one learns from this example is that imposing inequalities (45), (46) over the 
whole pp resonance range provides an approximate sufficient condition for a suitable field 
distribution. Minor deviations may be allowed for some x pp , since these inequalities were 
derived from the simplified formula (37) and a fixed neutrino energy. In any case, though, 
the main characteristic of a field increasing by a factor > 6 over a fraction shorter than 
9% of the solar radius seems essential, in order for the field to be compatible with solar 
neutrino data: the smearing effect due to the energy distribution is not as large with pp 
as with 8 B neutrinos, so this criterion should be retained. 

Let us reverse the attitude and search for a field distribution satisfying (45), (46) 
behaving linearly with x for simplicity. The 7 Be resonance is assumed to be located at 
the bottom of the convective zone where the field is taken to be maximal (x Be = 0.71, 
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A 2 m 2i = 1.6 x l(T 8 eV 2 ). It then falls to 2 x 10 3 G at the surface of the Sun. Taking 
/ = 10~ 12 , one has B(x7 Be )>6.6 x 10 5 G and (46) becomes 

B(x pp ) < 6.11 x itf—^G (59) 

with Epp given in eV. Given the fact that pp resonances start at 0.065 below X7 Be , one has 
essentially two kinds of options: 

1. First option. One imposes inequality (46) to hold only for part of pp resonances 
(x pp <X7 Be — Ax) thus allowing for its violation in the interval xt Bs — Ax<x pp <X7 Be — 
0.065. Consequently B{x) has to further decrease inwards from its value at x pp = x? Be — 
Ax in order to compensate for excessive neutrino conversion in this interval. (Recall 
that pp neutrinos are practically unconverted). So, taking for instance Ax = 0.09, the 
corresponding energy E pp in (46) for A 2 m2i = 1.6 x 10 _8 eV 2 is 3.17 x 10 5 eV. This gives 
B(x pp )<1.08 x 10 5 G from (59) {B(x pp )<9A x 10 4 G if one required all pp neutrinos to 
satisfy (46), i. e. Ax = 0.065). Hence, choosing 5(x pp )<9.4 x 10 4 G to allow for an 
overlap, the field will be 

B(x) = x<0.606 

B{x) = 6.73 x 10 6 (x- 0.62) + 9.4 x 10 4 0.606<x<0.71 (60) 

B{x) = -2.407 x 10 6 (x - 0.71) + 7 x 10 5 0.71<x<l. 

However the probability constraints on P B , Pi, Pl give respectively 

8 B (4.8 - 6.8) x 10- 9 eV 2 (1.95 - 2.70) x lO~ 7 eV 2 

7 Be (1.12 x 10~ 9 - 2.1 x 10~ 8 )eF 2 (61) 

pp < 8.5 x 10- n eV a 1.1 x 10~ 8 eV 2 , 

with fi u = 10~ 12 {i B , showing an inconsistency with 8 B. This is solved by requiring a more 
moderate slope of the field along the outer layers of the convective zone. A decrease from 
2.2 x 10 5 G to 2 x 10 4 G at x — 0.91 is sufficient while keeping the same values at the 
surface and X7 Be = 0.71. Thus the same suppression of 8 B neutrinos is achieved at larger 
densities than those corresponding to the interval (4.8 — 6.8) x 10 _9 ey 2 in (61). The field 
is then (see fig. 3): 



B(x) = x<0.606 

B(x) = 6.73 x 10 6 (x- 0.62) + 9.4 x 10 4 0.606<x<0.71 

B(x) = -3.4 x 10 6 (x - 0.71) + 7 x 10 5 0.71<x<0.91 

B(x) = -2 x 10 5 (x - 0.91) + 2 x 10 4 0.91<x<l 



(62) 



with a solution (/ = 10 



-12\ 



8 B (1.37 - 1.90) x 10" 8 e1/ 2 (1.95 - 2.70) x !0~ 7 eV 2 

7 Be (3.2 x 10- 9 - 2.1 x \Q- 8 )eV 2 (63) 

pp < 1.34 x" 10 eV 2 > 1.65 x 10~ s eV 2 . 

Hence A 2 m 2l = (1.65 - 1.90) x 10~ 8 eV 2 . 

2 It is worth noting that the range of possible neutrino magnetic moments extends in this case down 
to fx v — 6.5 x 10 _13 ^s- 
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2. Second option. Inequality (46) is imposed as from X7 Be — 0.065 inwards. B(x) may 
then ramain constant below this value since it will be sufficiently small throughout all 
x pp to ensure enough neutrino survival. Again a large slope at large x (upper convective 
zone) must be avoided for consistency with 8 B neutrinos. So B(x pp ) = 9.4 x 10 4 G for all 



pp 



and, as before, B(x7 Be ) = 7 x 10 5 G, (see (fig. 4)): 

B(x) = 9.4 x 10 4 x<0.645 

B(x) = 9.32 x 10 6 (a; - 0.645) + 9.4 x 10 4 0.645<a;<0.71 

B(x) = -3.4 x 10 6 (x-0.71) + 7x 10 5 0.71<x<0.91 
B{x) = -2 x 10 5 (x- 0.91) + 2 x 10 4 0.91<x<l 



(64) 



with a solution (/ = 10 



-12\ 



8 B (1.37- 1.91) x 10~ 8 eV 2 (1.73 - 2.37) x 10- 8 eV 2 

7 Be (3.2 x 10~ 9 - 1.95 x 10- 8 )eV 2 (65) 

pp < 1.34 x 10- 10 eV 2 > 1.5 x 10- 8 e1/ 2 . 

The range of solutions is slightly enlarged here relative to option 1: A 2 m2i = (1.50 — 
1.91) x lO~ 8 eV 2 . 

Such a large field at the bottom of the convective zone (7 x 10 5 G) is a consequence of 
requiring a comparatively small value of the neutrino magnetic moment (fi u = lO~ l2 fi B ) 
and, especially, the 7 Be resonance to to lie at 0.71Rs instead of lower densities as in the 
previous example (50). For lower densities, adiabaticity is in fact satisfied with lower 
magnetic fields (see eq. (25)). 

At this stage it should have become obvious that a monotonically decreasing field from 
the core or the inner radiative zone to the solar surface such as the ones proposed in refs. 



44], [42[ must be ruled out in the light of solar neutrino data. Also linearly varying profiles, 
although not to be excluded, can only provide borderline solutions. 

To conclude this section, four more field profiles rising sharply at the bottom of the 
convective zone will be briefly reviewed. The first two are variants of one of the fields 
proposed in ref. |43| namely, 



B(x)= ,^ M8Xl °\ 1 G , 0.7KK1 (66) 
v ; ch[20(x - 0.71)] ' v ; 

with B = 2.8 x 10 3 G and B = for x < 0.71 (see fig. 5). For the first of these one obtains 
the following two limiting cases 

8 B (1.60 - 2.50) x lO-W 2 (1.10 - 1.54) x 10- 7 eV 2 

7 Be (5.2 x 10~ 9 - 1.6 x 10" 8 )e1/ 2 (67) 

pp < 7.4 x 10- n eV 2 > 1.43 x 10~ 8 e1/ 2 

with \i v = 3.4 x 10~ n fi B and 

8 B (1.50 - 2.40) x 10~ 8 eV 2 (1.10 - 1.55) x 10- 7 eV 2 

7 Be (4.9 x 10~ 9 - 1.6 x \Q- 8 )eV 2 (68) 

pp < 6.9 x lO^eF 2 > 1.60 x 10~ 8 e1/ 2 

with fijy = 3.65 x 10~ n /iB- This leads to the range of solutions: 

3.4 x 10~ n /ii?<^<3.65 x IO^Vb , 1.55 x 10- 8 e1/ 2 <A 2 m 2 i<1.60 x lO^eV^ 2 (69) 
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with (46) satisfied for x vv <xi B e — 0.082. As opposed to the previous examples, the field is 
constant below its maximum, so the fraction of pp resonances that do not satisfy (46) has 
to be smaller. For the second of distributions (66), with the field limited to the convection 
zone, the range of solutions necessarily extends. In fact the two limiting cases are now 

8 B (1.60 - 2.50) x 10" 8 e\/ 2 (1.0 - 1.42) x 10~ 7 eV 2 

7 Be (5.1 x 1(T 9 - 1.60 x W 8 )eV 2 (70) 

pp < 6.4 x 10~ n eV 2 > 7.3 x !0~ 9 eV 2 

with fi u = 3.4 x 10 _11 /iB and 

8 B (3.3 - 7.3) x 10- 9 eV 2 (1.10 - 1.44) x 10~ 7 eV 2 

7 Be (2.0 x 1(T 9 - 1.6 x Kr 8 )e\/ 2 (71) 

pp no solution > 7.3 x 10~ 9 eV 2 

with n v — 1.3 x 10 -1 V B . The range of solutions is now 

3.4 x 10 -1 Vb<^<1-3 x IO^Vb , 7.3 x 10~ 9 eV 2 <A 2 m 21 <1.60 x Kr 8 e\/ 2 . (72) 

The last two profiles to be examined bear no discontinuity. They are 

B(x) = 2.16 x 10 3 G aKO.7105 

B(x) = 8.7 x 10 4 [1 - (2^5) 2 ]G 0.7105<x<0.7483 " ' ' 



with 
and 



B(x) = 10 & [1 - 3.4412(x - 0.71)]G 0.7483<x<l (74) 
8 684 x 1 D 4 

B(x) = , r ^ — rrG 0.7483<x<1 (75) 

v ; ch[20{x - 0.7483)] y ' 

shown in figs. 6 and 7 respectively. The first of these has a very limited range of solutions 
of which some violate up to 30% the inequality (45). Also (46) is satisfied only for 
E < 0.381MeV , a similar situation to that encountered with the linear field distribution 
(50). Again, this is a borderline case and the limiting solutions are 

8 B (6.2-9.1) x lO~ 9 eV 2 (8.6 x 1(T 8 - 1.2 x lO~ 7 )eV 2 

7 Be (3.5 - 6.5) x l(T 9 e\/ 2 (76) 

pp < 1.0 x 10- 10 eV 2 > 6.5 x l(T 9 e\/ 2 

with \i v = 4.1 x 10~ 12 fi B and 

8 B (4.8 - 6.8) x 10- g eV 2 (9.5 x 1(T 8 - 1.3 x lO~ 7 )eV 2 

7 Be (1.3 x 1(T 9 - 1.25 x l(T 8 )e\/ 2 (77) 

pp < 9.2 x 10- n e\/ 2 > 6.9 x l(T 9 e\/ 2 

with \x v = 6.15 x lO~ 12 fi B . These lead to the following range of solutions 

4.1 x 10 _12 /x B <^<6.15 x lQ- 12 fi B , 6.5 x 10" 9 e\/ 2 <A 2 m 2 i<6.9 x lQ- g eV 2 . (78) 
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In contrast to (73), (74), the second combination (73), (75), provides a much larger set of 
solutions whose limiting cases are 

8 B (1.29 - 2.1) x 10~ 8 eV 2 (9.5 x 1(T 8 - 1.3 x 10~ 7 )eV 2 

7 Be (4.33 x 10" 9 - 1.32 x l(T 8 )e1/ 2 (79) 

pp < 5.6 x 10- u eV 2 > 6.8 x lO~ 9 eV 2 

for u u = 7.2 x 10~ 12 fi B and 

8 B (3.30 - 7.32) x 10~ 9 eV 2 (1.0 - 1.4) x 10~ 7 eV 2 

7 Be (1.85 x 10~ 9 - 1.51 x 10- 8 )eV 2 (80) 
no solution > 7.32 x 10~ 9 eV 2 

for [l v = 2.1 x 10~ u fi B , shown in fig. 8. Hence, for the field (73), (75), the set of solutions 
is 

7.2 x 10" 1 Vb<^<2.1 x IO^Vb , 7.3 x 10" 9 eV 2 <A 2 m 2 i<1.3 x 10" 8 e1/ 2 . (81) 

The distinguishing feature between (74) and (75) is the same as between field profiles 
(60), (62): two fall offs along the convective zone, one linear and the other with an up- 
ward facing concavity. The second type is clearly favoured while the first is disfavoured. 
A downward facing concavity along the convective zone thus seems to be excluded. The 
type of exponential fall off proposed years ago [12] where this downward concavity pre- 
dominates, namely (see fig. 8), 

B(x) = — B \_ om , 0.7483<x<l (82) 
1 + exp^^ 

is also not possible. The main reason is that the moderate suppression required for 
8 B neutrinos demands their resonances to be mainly located in a zone where the field is 
substantially smaller than at the 7 Be resonance. But then the pp resonances are 'dragged' 
into a zone where the field is far too large. 

The main requirements for a suitable field, as from the analysis of this section are 
therefore: 

1) A very sharp rise across the bottom of the convective zone by at least a factor of 
6-7 over a distance not greater than approximately 9% of the solar radius. This sharp 
increase may be a discontinuity. 

2) A decrease along the convective zone whose slope is at first strong and then reduces 
to become very moderate or even closer to constant as one approaches the surface. 



5. Outlook and Conclusions 

The solar neutrino data were used to obtain information or, possibly, 'derive' the likely 
profile of the solar magnetic field on the basis of the existence of a neutrino magnetic 
moment and the RSFP mechanism. The main point of the analysis relies on associating 
neutrino energy with resonance location within the Sun. The experimental data and 
general basic features of solar models overwhelmingly favour an energy dependent neutrino 
suppression, the main uncertainty left being the amount of suppression of the highest 
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energetic S B neutrinos. So neutrino suppression is likely to be correlated to resonance 
location within the Sun and thus to the profile of the field intensity, since a strong field at 
the resonance is responsible for a high transition rate. These are the underlying arguments 
for a possible tomography of the solar magnetic field. 

On the basis of these arguments, the large suppression of the intermediate energy 
neutrinos consisting mainly of the 7 Be line at E = 0.86MeV, along with the total or 
nearly total survival of the detected pp neutrinos (0.236MeV r < E pp < 0.42MeV) is a clear 
indication of a solar field rising sharply along a relatively short distance. Such a topology 
can only be expected to be located around the bottom of the convective zone within a 
layer of at most (9 - 10)% of the solar radius. From the analysis of the previous section the 
field should rise by a factor of at least 6-7. All is consistent with this rise being extremely 
sudden, even in the form of a discontinuity, or to span most of that layer. If one believes 
in the RSFP mechanism, one may get an upper bound on the field (see eq.(46)) from the 
inner edge of this zone down to approximately 0.59i?s to be B max ~(l.l — 1.3) x 10 4 G, 
taking \i v = 10~ n fi B . At the peak the field should be B min ~7 x 10 4 G for the same \i v . 
Both these bounds are inversely proportional to the neutrino magnetic moment. Deeper 
into the radiation zone the field is essentially unknown, although there is no ground to 
believe in other large spatial variations. 

As regards the convective zone, one expects a decreasing field with no sudden vari- 
ations, the analysis favouring a larger slope towards the bottom with a more moderate 
one towards the surface. The 'true' field profile thus appears to be consistent with any of 
the ones shown in figs. 2 through 7, with fig.8 excluded because of the shape of the field 
decrease along the convective zone. 

In the original papers [[| introducing the hypothetical anticorrelation of active neutrino 
flux with solar activity as evidenced through the sunspots, it was suggested that the 
sunspot activity effect could extend as deep as the bottom of the convective zone. This is 
however by no means clear, hence the extent of neutrino trajectory affected by this 11 year 
sunspot cycle is unknown. The evaluation of the adiabaticity parameter 7c (eq.(25)) for 
all magnetic field distributions considered shows that only the one in fig. 5 (eq. (66) with 
B — for x < 0.71) may generate moderately adiabatic transitions for the Super-Kamio- 
kande neutrinos with energies ranging from their 7MeV threshold up to 9MeV. Even this 
situation will only occur at the upper edge of the solution, namely \i v = 1.3 x 10~ w [j jB , 
A 2 m 2 i = 1.6 x 10~ s eV 2 , (eq. (72)), in clear conflict with the astrophysical bounds on \x v 



45] . For all other cases the Super-Kamio-kande neutrinos are expected to undergo non- 
adiabatic transitions in the Sun, so no anticorrelation with sunspot activity is expected, 
even if the effect extends deep down into the convective zone. For the Chlorine experiment 
on the other hand, a larger fraction of S B neutrinos may be anticorrelated, since the 
energy threshold is much lower, but the effect may be totally unclear due to the dominant 
influence of the more energetic neutrinos whose resonances are non-adiabatic. Therefore 
RSFP scenarios are fully consistent with the absence of anticorrelation in the experiments 
so far |Lf| . 

As far as the approach used throughout this paper is concerned, a few words must 
be added regarding the simplified treatment of the intermadiate energy neutrinos. An 
overall suppression factor of at least 96% was considered for all of them, with a fixed 
energy E = Ei Be = 0.86MeV. They however span an energy range up to 1.7 MeV with 
this 7 Be line accounting for the largest part (~ 71%). Nevertheless, a more elaborate 
procedure is hardly justified at the present stage of both data and theoretical models and 
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should not change significantly the results derived here. In fact the suppression peak will 
necessarily be still assumed to occur at the 0.86MeV 1 Be line, with the other 1 Be line 
at 0.38 MeV, immersed in the pp energy range, being naturally much less suppressed. 
The CNO continuum may also suffer this suppression peak around 0.86 MeV but, for 
the remainder, it need not, along with the pep line (E = 1.442MeV), undergo a large 
suppression. 

It should be stressed that all information that can be obtained from the RSFP on 
the solar magnetic field only concerns its profile and not its absolute magnitude because 
the order parameter of the analysis is the product fiB. The analysis and results refer 
only to the transverse component of the magnetic field along the direction of propagation 
and nothing can be inferred as to its latitude dependence. They necessarily bear the 
uncertainties common to all solar physics models. 

From the experimental side, more abundant statistics and precise data together with 
new solar neutrino experiments are essential in order to improve our knowledge, along with 
independent laboratory searches for neutrino magnetic moment effects. In has been shown 
in this paper that the necessary fi v to provide a solution to the solar neutrino problem, 
thereby leading to information on the magnetic field, is within most of the relatively 
stringent astrophysical bounds |45| . Although this lies almost 6 orders of magnitude above 
the electroweak standard model prediction, one should keep in mind that the standard 
model itself is full of fine tuning and similarly large hierarchy problems affecting for 
instance fermion masses. Furthermore, there exist successful models |46] for large \i v 
which do not conflict with the necessary smallness of neutrino mass. 
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Experiment Ref. 



Data 



Units 



Homestake 
Kamiokande 
SAGE 
Gallex 
SuperKamiokande 







B 




2.54 ±0.14 ±0.14 
2.80 ±0.19 ±0.33 

76.2 ±6.5 ±5 

^• dl=c 0.13 



SNU 
10 fl cm -2 s -1 

SNU 

SNU 
lO^m" 2 ^ 1 



Table I. Neutrino event rates measured by solar neutrino experiments. 



TCL[[T6[ TCCCD[[T7[ P94[[Tg] RVCD96[[19| FRANEC96gg BP95[gT 



Chlorine 



pp 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


pep 


0.22 


0.21 


0.222 


0.221 


0.224 


0.224 


7 Be 


1.10 


0.995 


1.24 


1.15 


1.08 


1.24 


13 N 


0.063 


0.104 


0.109 


0.095 


0.0901 


0.105 


15Q 


0.21 


0.37 


0.379 


0.327 


0.306 


0.371 


8 B 


4.63 


4.06 


7.19 


7.03 


5.73 


7.35 


total 


6.23 


5.75 


9.1 


8.8 


7.4 


9.3 


Gallium 














PP 


71.1 


70.6 


69.7 


70.1 


70.7 


69.7 


pep 


2.99 


2.795 


2.99 


2.97 


3.01 


3.01 


7 Be 


30.9 


30.6 


37.9 


35.1 


32.9 


37.7 


13 N 


2.36 


3.87 


3.95 


3.45 


3.27 


3.82 


15Q 


0.21 


6.5 


6.46 


5.58 


5.22 


6.32 


8 B 


10.77 


9.31 


15.7 


15.4 


12.5 


16.1 


total 


122.5 


124 


137 


133 


128 


137 


Kam 
















4.4 


3.8 


6.48 


6.33 


5.16 


6.62 



Table II. Contributions of different sources of neutrinos to total capture rates in the Chlo- 
rine, Kamiokande and Gallium experiments from each standard solar model considered. 
All models include diffusion except for TCL [IE] and TCCCD [I7|. Units are 10 6 cm~ 2 s _1 



for Kamiokande and SNU for Chlorine and Gallium. 
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Fig.l Survival probability of pp neutrinos with a magnetic moment /i„ = 1.3 x 10~ 10 /i-B 
for the field distribution given by eq.(66) with B = in the radiative zone and core. The 
left end of the diagram corresponds to the resonance being close to the surface of the Sun 
or in the vacuum where the field is negligible, whereas the right end corresponds to a 
resonance in the interior (radiative zone and core) where the field also vanishes. 

Fig. 2 Field distribution given by eq.(50). Coordinate x denotes the fraction of the 
Sun's radius and the magnetic field is given in Gauss. The solution range is given by 
eq.(57). 

Fig. 3 Field distribution given by eq.(62) ('option 1'). For the solution see eq.(63). 

Fig. 4 Field distribution given by eq.(64) ('option 2'). For the solution see eq.(65). 

Fig. 5 Field distribution given by eq.(66) with B = below the convective zone. The 
range of solutions is greatly extended relative to previous cases (see eq.(72)) because the 
rise in the field at the bottom of the convective zone is both very large and sudden so 
that inequality (45) is largely satisfied. Also from the bottom of the convective zone to 
the solar surface the field has the 'adequate' shape that can account for a large range in 
terms of the 8 B solution. 

Fig. 6 Field distribution given by eqs.(73), (74). The range of solutions is very limited 
(see eq.(78)). 

Fig. 7 Field distribution given by eqs.(73), (75). Although the magnitude of the field 
is no larger than in the previous case, the solution is greatly extended (compare eq. (78) 
with eq.(81)) owing to the difference in shapes. 

Fig. 8 Field distribution given by eq.(82). Although this distribution has the correct 
shape across the bottom of the convective zone with a large and sharp increase, it fails to 
provide a solution, owing to its shape along the convective zone. 
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